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Abstract 

We propose an exactly solvable self-consistent kinetic 
model of polymerization-induced phase separation (PIPS) via spinodal de- 
composition. Using modified Cahn-Hilliard and Glotzer-Coniglio theories for 
early and late stages of spinodal decomposition, we find scaling regimes and 
compare the obtained results with existing experimental and theoretical in- 
formation on PIPS kinetics, finding a qualitative agreement. 
PACS: 05.70Ln, 61.41.+e, 64.75+g 
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I. INTRODUCTION 



In the last several years, polymer-dispersed liquid crystals (PDLCs) have generated sig- 
nificant interest among scientists because of their utility for display applications (see, e.g., 
[fj - [0]). Among the most important problems in PDLC preparation is to control the phase 
separation between the polymer and the liquid crystal. It is known that such phase sepa- 
ration can be achieved either by performing polymerization in the monomer - liquid crystal 
(LC) mixture at high temperature (where polymer and LC are miscible) and then quenching 
to low temperatures and thus inducing phase separation [temperature induced phase sep- 
aration (TIPS)], or by performing polymerization at such temperature at which monomer 
and LC are miscible, yet polymer and LC are immiscible [polymerization induced phase 
separation (PIPS)]. In the former case, the phase separation is controlled by the quench 
depth, in the latter case, it is controlled by the polymerization kinetics. 

Recently, several studies of PIPS kinetics, experimental ( || - flOU), theoretical [p| , and 
numerical ( [jTJ] - |13|), have been performed. These results showed that PIPS is reasonably 
well described by theories of ordinary phase separation induced by a thermal quench, so that 
at the early stage of phase separation (PS), classical linear theory of Cahn and Hilliard ( |14| 
- [|16|) is applicable for the initial growth of phase regions; at later stages, scaling models of 
Binder and Stauffer |1| , Langer, Bar-on and Miller [17 and Siggia []TU[ are used to describe 
the behavior of scattering function or domain size. At the same time, no model exists yet 
that describes the analytical dependence of PS kinetics on polymerization rate. The purpose 
of this work is to provide at least a qualitative approach to the problem. 

The paper is organized as follows: in section |I| we introduce thermodynamic and kinetic 
variables used in the problem, discuss the phase diagram and make approximations about 
conditions of polymerization. In section |T| we derive hydrodynamic equation of PS and 
discuss the self-consistent approximation. Then, in section [IV|, we analyze approximate 
scaling solutions of these kinetic equations and compare them with existing experimental and 
theoretical information. Finally, in section [V], we discuss the validity of the approximations 
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used and analyze possible practical use of the model and future changes and improvements. 



II. PHASE DIAGRAM AND THERMODYNAMIC VARIABLES 



We will consider a system consisting of monomer, polymer and LC in isotropic state. 



Following Teixeira and Mulder [E|, we will neglect polydispersity of polymer lengths and 
disregard any difference between gelation and polymerization. In this case, the system under 
consideration can be described as a ternary mixture with being the volume fraction of LC, 
M being the volume fraction of the monomer, and P being the volume fraction of polymer. 
We assume an incompressibility condition that reads: 

<p + M + P = l. (1) 

Because of condition ([[]), only two of the three concentrations are independent, and it is 
more convenient to work with variables <fi and Y = P/(P + M) - the degree of conversion. 
We assume that polymer and monomer are always miscible, and that the monomer - LC 
interaction and the polymer segment - LC interaction are the same. In this case, the Flory- 
Huggins free energy of the mixture can be written as: 

F = J dy[01n0+(l-r)(l-0)ln(l-0) + X 0(l-0) + ^(V0) 2 ], (2) 

where \ is Flory-Huggins parameter, related to temperature (it is usually assumed that 
X = A + B/T, where B is a positive constant, A is an arbitrary constant). 

The schematic phase diagram for the system (Q) is shown in Figure 1, with different 
curves showing spinodal lines for different values of conversion parameter Y . It is similar 



to spinodals of polymer-LC system derived in p0[ . It can be seen that with Y increasing, 
the critical temperature and critical concentration move toward high temperatures and low 
concentrations of polymer, respectively. For each value of the Flory-Huggins parameter 
X between Xci and Xc2, there exists a critical concentration <p c (x) an d a critical degree of 
conversion Y c (x), so that the point (x, 4> c (x)j ^c(x)) * s a critical point. 
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We will assume that polymer - LC phase separation occurs only via the spinodal decom- 
position (SD) mechanism and neglect nucleation in the metastable phase. Further, we will 
assume that the temperature at which the process occurs, is not very far below the critical 
polymer-LC temperature. In that case, as it will be shown later in this section, expression 
for the free energy used in deriving kinetic equations, can be simplified dramatically. 

For a given x, the critical composition (f> c (x) an d the critical conversion Y c (x) are: 



Y C ( X ) = V 8 * " 2 X, (3) 

«■ = /J- « 

The isothermal (fixed x) spinodal curve, such as shown in Figure 2, can be described by 
the following equation: 



, , Y-Y (Y -Y C )(Y + Y C + A X ) Y-Y c 
*-*•*>— ±\/ 16? < 5 > 

the second equality valid in the limit Y,Y C « 4%. We then construct the Ginzburg-Landau 
free energy density corresponding to this modified spinodal (i.e., all minima of this free 
energy expression lie on the spinodal): 

faL = - Y -^* 2 + XX\ (6) 
[3F = |^[/ GL + ^(V0) 2 ], (7) 

where X = <fi — <p c . 

In the next section, we will derive kinetic equations and analyze their solutions in the 
self-consistent approximation. 



III. KINETIC EQUATIONS 



In order to fully describe PIPS kinetics, it is necessary to determine time dependence of 
both degree of conversion Y and phase separation X. The former variable is a nonconserved, 
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the latter one is a conserved order parameter. Assuming that polymerization occurs as a 
simple bimolecular reaction with a rate /C, we can write equation for Y as follows: 

^ = /C(l-F) 2 (l-0 c -X) 2 . (8) 

For the phase separation, we employ the classical time-dependent Ginzburg-Landau 
(TDGL) equation: 

§ = rv^-KV^]. (9) 

As written, equations (H) - @ constitute a very complex system of coupled nonlinear 
partial differential equations, susceptible only to limited numerical treatment. We make 
several approximations in order to somewhat simplify them and make them suitable for an- 
alytical study. First, we average spatially the equation (§) to obtain a spatially independent 
Y{t) =< Y(r,t) >: 

dY 

^ = /C(l-F) 2 [(l-0 c ) 2 +<X 2 >]. (10) 

The system of coupled equations (|| - |10D somewhat resembles the kinetic equations that 
we derived for the annihilation-diffusion systems: 

f t = -/C(p 2 - < f >)), (11) 



dt 



DV 2 /(r,t)=0, (12) 



where the (averaged) particle density p is a nonconserved order parameter decaying to 
during the annihilation process, and the charge density / is a conserved order parameter 



whose fluctuations play important role in the annihilation process ||21|| . However, the PIPS 
kinetics is more complicated than the annihilation-diffusion process because equation (|9]) is 
nonlinear in X. On the other hand, it seems that in polymerization, fluctuations play a 
less important role than in annihilation, and it is possible to decouple the two by neglecting 
the < X 2 > term in (|i0|) and essentially having Y(t) a known function to be plugged in 



equation (§). Therefore, we will not consider equation fllOD in our subsequent analysis and 
instead concentrate only on the equation (|). 

We will also assume critical polymerization, i.e., that the system enters the two-phase 
region on a phase diagram through a critical point, undergoing a second-order transition. 
Mathematically, this means that < X(r,t) >= 0. If this is the case, PIPS occurs in two 
stages: first, when the concentration of polymer is small, and the mixture is still in the one- 
phase region, kinetics is diffusive, fluctuations are generated spontaneously due to random 
forces, and gradually shrink due to diffusion; second, when the degree of conversion becomes 
larger than the critical value Y c , kinetics changes from diffusion to spinodal decomposition, 
and fluctuations begin to grow rapidly. This (rather trivial) observation is supported by 
experimental evidence that PIPS has a relatively large induction time. 

In the remainder of this section, we will study in detail the phase separation stage of 
the process. Let us rewrite equation (|^) in Fourier space and simultaneously transform 
everything to dimensionless variables: 

= < x\t) > +e}x(Kt) (is) 



where X — ► XyJ^, k — ► k^j^r, t — ► t r(1 / e) . We also substituted X 3 with 
< X 2 > X, as is usually done in self-consistent theories ( P^fl , P3fl , P^fl ). Equation (|I3"D can 
be readily analyzed if one knows the t-dependence of Y. This analysis will be done in the 
next section. 



IV. SOLUTIONS OF KINETIC EQUATIONS - TWO REGIMES 

A. Early Stage of Spinodal Decomposition 

In the early stage of spinodal decomposition, fluctuations are relatively small, and it is 
possible to linearize the equation flTBD by neglecting the < X 2 > term. Such linearization is 
used in the original Cahn-Hilliard theory of phase separation after quench into the two-phase 
region. In PIPS case, however, the system moves into the two-phase region as polymerization 
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proceeds, and so the transition happens at some time to > 0. Here we will not attempt to 
determine t , but in principle, it is not difficult if polymerization rate is known: at t — t , 
Y = Y c . For the early stage of the process, then, we can write: 

\^ = JC(t-to). (14) 

At later times (t — t > /C _1 ), we just assume \Zy° = 1- 

The linearized TDGL equation has as its solutions all excitations of the form: 

X(k, t) = X(k, t ) exp[u k {t - t )], (15) 
with the dispersion relation: 

= 2 ^^|) " k ") = ^(^^ - k 2 ). (16) 
Excitations that grow most rapidly have a wavevector k m given by: 



k m = (17) 

Unlike the thermally induced phase separation, where k m is constant during the early 
stage, in PIPS it changes with time even in the early stage. The fact that k m increases with 
time indicates that the predominant mechanism is the formation of new fluctuations rather 
than the growth of existing ones. 

At time t — to — /C _1 , the system have polymerized completely and yet it is still in the 
early stage of SD. From now on, the regular framework of SD theories is applicable. In 
the next subsection, we will use the self-consistent approximation with the scaling ansatz to 
describe the late stage of the SD in PIPS. 



B. Late Stage of the Spinodal Decomposition 

In the late stage of SD, we can no longer neglect the < X 2 > term, and equations of 
evolution become strongly nonlinear. Even though the self-consistent approximation allows 
iterative treatment, its complicated character prevents us from writing an exact solution 
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in analytical form. Coniglio and Zannetti |22| described the solution of TDGL in this self- 
consistent approximation to find that k m oc (t — to) -1 ^ 4 ( we omit logarithmic correction 
to the power law). Such behavior is called multiscaling and is exact only for the infinite- 
component order parameter; for the scalar order parameter, one should expect a slightly 
different exponent. The other scaling theories of PS, based on different mechanisms of 
domain growth, proposed exponents ranging from 1/3 (Lifshitz-Slyozov law, based on a 
coagulation model) to 0.21 (Binder and Stauffer model Hl8|), so the mean-field value of 0.25 
seems to work rather well. We will outline another way of solving the TDGL in self-consistent 
approximation, and in the process will calculate the scaling function for the structure factor. 

Let us start from the scaling ansatz for the structure factor S(k,t) =< 
X(k,t)X(-k,t) >: 

S(k,t) = k- d F(k/k m ), (18) 

where T is a scaling function which has to be determined. Since 
<X\t) >= j*fasQt,t), 

equation ( |T3"D can be rewritten in terms of T as: 

- dk-J-^T - k-^x^ = -2A£V(-1+ <X 2 > +k 2 m x 2 )F, (19) 

where x = k/k m . The scaling function T is time independent, and from this requirement 
we determine the critical exponent A ik m oc it — t )~ x ). First, we separate variables in 
equation (|T9|) to obtain the implicit expression for T: 

= _ d + Ml( ^ ) -V(-i + < » + (20) 

From the requirement that function T has a maximum at x = 1 (stemming simply from 
the definition of k m ), we obtain: 

- d + 2kU^)-\-l+ < X 2 >) + 2kl(^r = 0. (21) 



Moreover, all three terms in equation ( |2"T| ) should be time-independent, otherwise the 
scaling ansatz is not applicable. This immediately dictates the critical exponent A equals 
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1/4, as it is required to be in self-consistent theories of TDGL (see, e.g., Coniglio et al. [p2fl). 
If we attempt to make k m continuous by connecting the early stage asymptotics and the late 
stage asymptotics at point t — to — K>~ 1 -, we obtain the zero-order approximation for k m (t) 
, < X 2 {t) >, and F(x): 

k ™ ~ lQJC(t-t o y (22) 

< - Y>Kl -vr^ ; (23) 

'""^ ~ + (24) 



with normalization: 



dlnx 2JC 1 2/C 



Equation (^4|) can be solved to yield the explicit dependence of T on x: 

T = ^- d exp[-^ + (d + (26) 

where Tq is simply a normalization constant. 

In Figure 3, we compare our scaling function (|26| ) with the one experimentally measured 
by Palffy-Muhoray et al. . It can be seen that the proposed scaling function describes the 
long-wavelength part of the spectrum reasonably well, while in the short-wavelength region, 
it decays much faster than the experimental scaling function. This behavior is somewhat 
expected because the mean-field Landau-Ginzburg treatment fares poorly when fluctuations 
are large and/or domain walls are too steep. An oversimplified treatment of polymerization 
process may also contribute to this discrepancy. However, agreement in the long- wavelength 
region, achieved with the help of just one adjustable parameter /C, is encouraging. 

Obviously, solutions (|22| ) - ( p^ ) can be considered only the zero-order approximation 
even in the self-consistent model. Nonlinear differential equations of evolution should be 
solved iteratively, and in this case, we expect to find logarithmic corrections to k m (t) and 



significant changes in the behavior of < X 2 >. Calculated self-consistent solutions thus 
should be viewed as a qualitative, rather than quantitative, description of the PS process 
and the influence of polymerization on its kinetics. 

The above-described late stage of the SD corresponds to a truly fractal system. It can 
be viewed as a set of LC-rich droplets in a polymer-rich matrix; within each droplet there 
are smaller polymer-rich droplets, and so on. This pattern also has been observed in the 
PIPS experimental studies [f|. 

At even larger times, the fractal picture is destroyed. As droplets grow larger, their 
interaction and overlap, which is not taken into account in the above treatment, becomes 
more and more important. The droplet size grows not only due to the SD mechanism, but 
also (and most importantly) due to aggregation and coagulation. Such a mechanism was 
described by Siggia |19j in his scaling theory, which yields A = 1 at this very late stage of 
phase separation (such behavior was indeed observed for PIPS in the same experiment ||.) 
However, this stage is not considered in our theory, which was primarily aimed at analyzing 
the influence of polymerization on the phase separation kinetics. It is rather obvious that 
at the time when Siggia mechanism becomes dominant, no polymerization takes place, and 
phase separation occurs in the same way for both PIPS and TIPS. 



V. CONCLUSIONS 



We analyze in a self-consistent approximation the critical polymerization-induced phase 
separation (PIPS). Contrary to the temperature- induced phase separation (TIPS), in the 
critical PIPS the phase separation starts only after some induction period, when polymer- 
ization advances enough to make liquid crystal strongly incompatible with the polymer. By 
assuming a stepwise polymerization rate and neglecting the feedback between phase separa- 
tion and polymerization, we simplified the problem and obtained analytical solution for the 
early (Cahn-Hilliard) and late (Coniglio-Zannetti) stage of the critical spinodal decomposi- 
tion (SD). 
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Among the other approximations: the isothermal process occurs at relatively high (close 
to the miscibility of polymer and liquid crystal) temperature - this approximation justifies 
the use of symmetric Ginzburg-Landau free energy; constant, rather than concentration- 
dependent, mobility parameter T; gelation is deemed unimportant or non-crosslinked poly- 
mer is considered. All these features make model very simple in comparison with real PDLC 
systems; at the same time, it obviously enables us to obtain a better insight into this com- 
plicated process using analytical, rather than purely numerical, method. 

Obvious goals for the improvement of the theory would be: expanding the model for the 
case of non-critical PIPS; allowing for the more realistic description of the polymerization 
process; accounting for the difference between gelation and polymerization. These studies 
will be undertaken in the future. 
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FIGURES 



Figure 1. Spinodals of polymer-monomer-liquid crystal mixtures. Each curve corre- 
sponds to a given degree of conversion from monomer to polymer, from Y = (no polymer) 
to Y = 1 (no monomer left). The solid horizontal line is a cross-section \ = 1-5. 

Figure 2. Spinodal points for polymer-monomer-liquid crystal mixtures at fixed temper- 
ature (x = 1.5). 

Figure 3. The scaling function F(q/q max ) vs the scaled wavevector q/q ma x- The data 
(six differentsets of points) are reprinted from Ref. |J; the solid line represents the best fit 
for these data using Eq. fl26|) with adjustable parameter /C = 0.025. 
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